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q ■ Abstract. The simplified models of interaction of charged matter with resonance modes of 

| radiation generalizing the well-known Jaynes-Cummings and Dicke models are considered. 

It is found that these new models are integrable for arbitrary numbers of dipole sources 
and resonance modes of the radiation field. The problem of explicit diagonalisation of 
corresponding Hamiltonians is discussed. 
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1 Introduction 

In the last four decades, much attention has been paid to the problem of describing processes 
of interaction of charged matter with electromagnetic radiation at resonance ^ |2J El El • This 
problem, being overcomplicated from the mathematical viewpoint in the general case, admits 
some very attractive simplifications which allow one to construct some rather simple models 
having even exact solutions. These models can be used for quantum statistical description of 
real interaction processes. 

The starting point for constructing such models is the quantum Hamiltonian 

^ ' H = ^-(p- e -A) + u(x) + Hp, Hp = £> fe (b+b k + , 



;_i ■ where p and m are the momentum and mass of the electron in some atom, u(x) is the Coulomb 

potential of interaction of the electron with nucleus, and A is the vector potential of the secondly 
quantized electromagnetic field which is given by the formula 

1/2 



where are polarization vectors and b^, b^ are the usual operators of annihilation and creation 
of the mode with wave vector k and the energy ojk- There are several steps of simplification of 
this general Hamiltonian. First, one neglects the term which is nonlinear in the vector potential, 
2^j^4 2 - Second, one considers only transitions between higher levels of an atom, 

p 2 + u(x) — ► u>oa z + const, 

2m 
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where a z is the Pauli matrix. Third, one considers only small k, ka <C 1, where a is the 
typical size of an atom, e ±lfcr ~ 1. Fourth, one uses the dipole and so-called rotating wave 
approximations [3], which result in 

e 

mc 



where <r ± are the usual Pauli matrices. And finally, one considers only the case of resonance: 
Sfc — * ^fixcd- By all these steps one can write down very simple Jaynes-Cummings Hamilto- 
nian [2] which is linear in the operators of creation and annihilation of one resonant mode, 

Hjc = oo b + b+^uj a z +g(a + b + a~b + ). (1) 

The solution of eigenequation for the operator is almost trivial, but it allows one to de- 
scribetime evolution of the corresponding wave function in great detail, with transparent ap- 
plications to physics of one-mode laser. This extremely simple model admits almost evident 
generalisations in two directions. 

2 The simplest generalisations of the JC model: 
the n-level atom and Dicke model 

The first of these generalisations is the model of one n-level atom with (n — 1) modes of resonant 
radiation, 

n—1 n—1 

Hjc = ^2 UJa ^ a + AR nn + ^2 g a (b a R na + b^R an ) , (2) 

ot=l a=l 

where 

N a = b^b a — Rata, {Rap) 1& = 8aj$/3S-> 

the sets {to a }, {g a } are energies of the modes and their coupling constants, and A is so-called 
detuning parameter |H] • The extended physical motivation for this kind of the generalization 
of the Jaynes-Cummings model can be found in There are evident commutation relations 

[N a ,N p ]=0, [H^,N a ]=0, 

which make the problem of diagonalisation of @ at n > 3 quite easy. The second way of 
generalisation of the Jaynes-Cummings model consists in considering an arbitrary number N of 
two-level sources of radiation, 

N 

H D = ub + b + + K a + b + v ( -b + )] ■ (3) 

3=1 

The solution of this so-called Dicke model which has been proposed many years ago is 
far from being trivial. It turned out that the problem with the Hamiltonian @ has some 
hidden dynamical symmetry which allows one to construct some nontrivial operators commuting 
with ©. 

What happens if one will try to use both the ways of generalisation of the Jaynes-Cummings 
model, i.e. consider the Hamiltonian of the problem of an arbitrary number N of multilevel 
sources interacting with an arbitrary number n — 1 of resonant modes of radiation, 

n-l N Ti-l N 

H ^ = y: u a N a + £ A k m +ye 9 i k) [b*m + , (4) 

a=l k=l a=l k=l 
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where N > 2, n > 3 (the detuning parameters {A} and coupling constants {g a } are assumed 
to be arbitrary)? The answer is not clear till now. In the following two sections we will argue 
that there is some simplification of the spectral problem defined by (0} in the case of equality of 

all the modules of the coupling constants {g}: these constants should not depend on k and their 

(k) 

dependence on a is given by the formula g a oc e a , s a = dbl. Namely, we will construct in this 
case explicitly some set of N — 1 nontrivial operators commuting with (jl} ; their existence allows 
one to call the model @ quantum integrable under the above choice of the coupling constants. 



3 The models of Gaudin type 

In this section, let us consider some more general operators introduced first by Gaudin [7j, 

n N+l 

n f ] =EE fjfiiffia* j = l,..., N + l, (5) 

a,/3=l sjtj 

where |^ Q/ g| are operators of absolutely arbitrary representations of SL(n), obeying the com- 
mutation relations 

[Ti^K]= 5 3s{nM- f U 5 ^)^ ( 6 ) 

and {fjf} are arbitrary sets of numbers, obeying the evident restriction f^f = /J^\ It is easy 
to see that operators 

N+l 

N„ = f j 

3=1 

trivially commute with 

[nf\N a ] =o. 

Now let us pose the problem: how to find such the sets of numbers {/} that 

[nf\n { k N) ]=o (?) 

for all indices {j, k}. It is easy to see that in this case the "Hamiltonian" 

N+l n 
j=l a=l 

where rjj and io a i v ha are arbitrary constants, commutes with all 7~L^ ', 

[Hi N \n ( j N) ]=o. 

In the simplest nontrivial case of n = 2, one can always choose {/} such that = —f^, and 
the equations ® can be written in the form 

N+l n 
sj^j,k o,/3,A=l 

x (T^xA ~ TijLTlp) = 0, (9) 
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which immediately gives a overdetermined set of bilinear equations for {/}, 

fjff£ + + fgfg = 0. (10) 

It has been shown by Gaudin |Sj that for n = 2 there is some so-called anisotropic solution 
to JTUJ), 

where are arbitrary numbers. Now, following Gaudin, let us choose them as 

7T 

Mfc = -eAfc + - +tijv+i, fc = l, ...,N. 

As e — ► 0, then ff% + i ~ eAj, f*N+i ~ 1- Let us cnoose {^g +1 } as Jordan-Schwinger realization 
of SL(2), = —b a bp and consider the limit A^ 2 3> 1, where A2 is an eigenvalue of the 

operator N2. It is easy to see that in this limit 6-operators can be treated as c- numbers; 

-1 

62, b~2 ~ \Ni\ let us also take e ~ N 2 2 and send N2 — > 00. In this limiting process, one easily 
finds 

N 

lim H™ = u lb th + [A, (fff - f + (ifffe + f W^)] • (12) 

iV2 — >00 L — * 

i=i 

And finally, if we will choose {2^jg } as cr-representation of SL(2), we will obtain just the Dicke 
Hamiltonian! 

This derivation of the Dicke model proposed first by Gaudin [Jj is not original of course. There 
are some recent papers in which the connection between it and the Gaudin models was exploited 
for construction the algebraic Bethe ansatz method for eigenvectors of the Dicke Hamiltonian 
and some generalisations which might be considered as the interaction term without the rotating 
wave approximation pi lTDlITT] . 

The above limiting procedure, however, does not work in the general case n > 3. One also 
has in this case the overdetermined system of bilinear equations 

fa/3 raX _i_ f a\ f/3\ , fCt/3 r-. 

Jjs J sk "T Jkj J js "r" J jk -Iks ~ U ' 

but, if n > 3, there is a possibility of the choice of indices a / ^ / A which excludes all 
"anisotropic" solutions of the type (jlll) . There are only "isotropic" solutions of the type f"" = 

= 7, -,, ; an d hence there are no parameters for N n — > 00. Hence it seems that the Gaudin 
approach does not lead to any integrable system of the type (j3J). 

Let us prove absence of "anisotropic" solutions to (|1()[) as iV = 2, n > 3. In this simplest 
nontrivial case the set {/} can be parametrized as 

333 
fjs = ' S ^ l £ ispA p P , f sk =^2 £ skpA p , fj 8 =^2 £ jspAp . (13) 
p=i p=i p=i 

As j ^ k ^ s, HOI) reads: 

AfAf + AfA p k x + AfAf = 0, 

with general solution 

AaX _ A . aX _ A A A 

smz/i smv2 sm(y\-\-V2) 
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A^ a = A cot vx , A% a =A cot v 2 , A% a = -A cot (i^ + u 2 ), 
where A, v\, v 2 are some parameters. Now, if one takes a 7^ /3 7^ A, one gets from (jlUj) 
.4 2 



2 cos » cos f cos 



,4 = 0. 



Surprisingly enough, one can construct integrable models of Gaudin type at n > 3 by using only 
"isotropic" solution of (|10|). It will be done in the next section. 

4 New integrable model with AT sources and n — 1 modes 

The receipt of constructing the integrable model of the type (jlj is quite unusual. Let us choose 
the above W-operators as obeying the equality 

[nf\ni N) ]=s jk ({f^},{x}), (i4) 

where x is some parameter and Sjk have the following property: as x — > 00, there are only zero 
eigenvalues of them in a subspace in which {rfgj act nontrivially. Then we can introduce new 

operators ( lim TiS ) in this subspace which will commute, and H^ +1 = — Y] TiS is some 
nontrivial Hamiltonian. By direct calculation of the commutator of the operators ©, one finds 

N n 
s^j,k a,f3,X=l 

n 

v (rfj T'k rjlS rpj rpk rpS \ ST^ ( f(X/3 raX j_ fCtX f/3X , fOc/3 fj3\ \ 

X K 1 a/3 1 Xa 1 /3X ~ 1 /3X 1 aX 1 X/3) ~ 2-> \Jj,N+lJN+l,k + Jkj Jj,N+l + Jjk h,N+l) 

a,P,X=l 

x (fi?TLT» +1 - q&it? 1 ) ■ (15) 

Let us choose /j^ = A(i/j — z^) -1 , 1 < j,k < N, i.e. use the isotropic Gaudin solution at 
j,s,k < N. Then the double sum vanishes and we are left with 

n 

Sjk{{Taf3}, {x}) = ~ ^2 [f^N+lfN+^k + _ {fk,N+l ~ fj,N+l)] 
a,(3,\=l i h 

x (TipTxJ^ 1 ~ T^ x f^) . (16) 
There are no other restrictions to the numbers except symmetry. Hence we have 



free parameters at given {^1, {z^}} in the relation Q16JI . Of course, there is trivial case 



Nn(n+1) 
2 

if all /"n+i = 0- Fortunately, one can construct the nontrivial ansatz: 

f^N+l = (PpnffjJ+l + Sanfj^+i) (1 ~ S an 8p n ) /jjv+1' 

i.e. fjpf + i = unless a ^ n or (3 7^ n. 



The commutator [Hj 1S s ^ too complicated. Let us make the next assumption: 

f$r +1 = e x VAF, f^ +1 = -Fuj, \e x \ = 1, (17) 
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where {vj~\ are arbitrary numbers. Then ()16|) can be recast in the form 



n-1 



/3,A=1 
n-1 



A=l 



Fvuf^(T j T k T N+l -T J T k T N+1 \\ 

r VkE\{± nX -L nn -L Xn ± Xn ± nn ± nX )]. 



(18) 



Let now choose 



T 



N+l 



as Jordan-Schwinger representation of SL(n)). Hence the basis of the whole Hilbert space of the 
problem consists of the direct products of the basis vectors of the spaces in which Tl x , 1 < j < N 

act, and various basic vectors of above Jordan-Schwinger representation, { {/} (^/?) ^|0)}> 
{1} £ N. Consider the action of the right-hand side of (17) on the subspace spanned by the 
vectors t/4m with l n = L, L S> 1 such that 



N m( L ) ~ r, 2 ( n( i ) /v <A L ^ 



N 



Then 



1 / j nr, 



bnK 



N 



-L 



f N+1 ~ 



1 < a < n - 1. 



If 



V- AA fnn f 

r ~ ~Y?' Jk,N+l — JY~ ' Jk,N+l 



An 



L 



then 



lim 

L— >oo . 



W = 



and 



wf = lim 

J L— >oo J 

commute on this subspace! The explicit form of 7^-operators reads 

N n 

1 ^ (,.. „ \-l , „^ 



s^j a,P=l 
n-1 

E^(^A + ^ n )+A^ 

A=l 

(JV) TV Tn-1 

i=i j=i fe=i 



(19) 



2 ex^{bxT* x + 6+f A fc J + A^fl 



A=l 



A? 



They commute with [ii, - } and ./Vq, = — Tqq, + 6+6 
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And finally, let us choose as T 3 „ the matrix representation of SL(n), 



f 3 -> B J 



and add linear combination w a JV a to H N ^ V 

a=l 

It is easy to see that we obtain just the operator of generalized Dicke model (JIJ) for n > 3, 

n-l JV n-l AT 



h <*»> = e u a N a + e + EE ^ + 6. 



a=l k=l a=l k=l 

with coupling constants 

gW = -Ae a VA, e a = ±1. (20) 

Hence, under the above choice of the coupling constants, the most general model (jlj becomes 
quantum integrable and the spectral problem might be simplified. However, we did not find 
the way for doing it except the simplest nontrivial case of N = 2 which is described in the next 
section. 

5 Explicit solution for N = 2 

In this case the Hamiltonian (@J) under the conditions l|20j) can be written as 

n-l n-l 2 

= £ Wa iV a + A Xj RL + A 2 R 2 nn + E E ML + b « k *n] > (21) 

a=l a=l fc=l 

2 

iv a = &+6 a -E^L- 

i=i 

Let us rewrite (|21(1 as 

n-l 

H n ] = ^2^ a N a + k 1 + k 2 , (22) 
a=l 

n— 1 n 

k,2 = ± (Ai - A 2 )" X A,2 + E Mna + + ^1,2^, A,2 = E ^ A- 

a=l a,/3=l 

Eigenfunctions of 1)22(1 are of course common eigenfunctions of (n + 1) operators {iV a }, fci, foj 

n— 1 n— 1 



A,u=l A=l 



Z„|0>, 



where <^ are eigenvectors of {j? AA }, {RiA and Z n = \\ (6+) a . 



MM J 

a=l 



One gets from the relations kiijj = kitp, k 2 ip = k 2 ip the following set of algebraic equations 
for the coefficients {^4, B, C}, 

A uA (A 1 -A 2 )- 1 + 4 2 ) = A; 1 A AAt , 
(A l -A 2 y 1 + C = hB^\ 
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B^ ) (A 1 - As)" 1 + A^f } +J2A flX (N fl + 1 + J AA4 ) = fcjjsf , 
C(Ai - As)" 1 + CAi + B ^(N„ + 1) = fciC, 



-A lxX (A 1 -A 2 )- 1 + B^ = k 2 A XfJL , 
-B { V(A 1 -A 2 r 1 + C = k 2 Bf\ 



-4 2) (A! - As)" 1 + A 2 B^ } +J2A Xlx (N fl + 1 + 5 Xfl ) = k 2 B x 1} , 

n 

-C(Ax - As)- 1 + A 2 C + Y1 B ? ] + X ) = 
They can be easily reduced to two equations for fci , A;s : 

n-1 

(fci + k 2 - Ai)(Ai - As)- 1 + (Ax - As)" 2 - h(h - Ax) + 1 + ^(iV^ + 1) = 0, 

n-1 

(A 2 - h - fe)(A 1 - As)" 1 + (Ax - A 2 )- 2 - fc 2 (A; 2 - A 2 ) + 1 + ^(iV M + 1) = 0, 
and the eigenvalues of the Hamiltonian (|21|) are given by 

n-1 



/i 2,n = ^2uj a N a + k! + k 2 . 



a=l 

6 Conclusions 

In this paper, we found the most general quantum integrable model of interaction of (n — 1) 
modes of radiation with N dipole sources which comprises all known models of that type: Jaynes- 
Cummings model (N = 1, n = 2), its generalisation for arbitrary number of modes (N = 1), 
Dicke model (arbitrary N, n = 2). The commutative ring of operators which includes 
is found in explicit form. We did not use any "anisotropic" form of the Gaudin solution; the 
Hamiltonian of the model was obtained via some limiting procedure. We almost immediately 
got solution for N = 2, arbitrary n with the use of additional integrals of motion, but the 
most interesting case of solution for arbitrary N and n is unreachable at the present stage of 
investigation. We hope to come back to this problem in the future. 



This paper has been presented for memorial volume of Vadim Kuznetsov. We both knew him 
personally — he was our guest in Dubna almost 17 years ago when he was a PhD student of 
Professor I. V. Komarov in Leningrad. We remember him as bright young man, very active in 
the field of quantum and classical integrable systems. Later he solved very complicated problem 
of integration of equations of motion for classical Toda chains with non- exponential ends |12j 
proposed by one of us (V.I.). 
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